Quantum Hall Effect in Biased Bilayer Graphene 
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We numerically study the quantum Hall eflFect in biased bilayer graphene based on a tight-binding 
model in the presence of disorder. Integer quantum Hall plateaus with quantized conductivity 
(Jxy = ve^ jh (where v is any integer) are observed around the band center due to the split of the 
valley degeneracy by an opposite voltage bias added to the two layers. The central (n = 0) Dirac 
Landau level is also split, which leads to a pronounced = plateau. This is consistent with the 
opening of a sizable gap between the valence and conduction bands. The exact spectrum in an 
open system further reveals that there are no conducting edge states near zero energy, indicating an 
insulator state with zero conductance. Consequently, the resistivity should diverge at Dirac point. 
Interestingly, the = insulating state can be destroyed by disorder scattering with intermediate 
strength, where a metallic region is observed near zero energy. In the strong disorder regime, the 
Hall plateaus with nonzero v are destroyed due to the float-up of extended levels toward the band 
center and higher plateaus disappear first. 

PACS numbers: 73.43.Cd; 73.40.Hm; 72.10.-d; 72.15.Rn 



I. INTRODUCTION 

The discovery of an unusual quantum Hall effect 
(QHE) in bilayer graphene has stimulated great inter- 
est in the study of the electronic transport properties of 
this new material [E i, i, i, d, @, 0, i, %, E, Ell, El, [H ■ 
At low energies and long wavelengths, the electrons in bi- 
layer graphene can be described in terms of massive, chi- 
ral, Dirac particles. While previous studies have focused 
on unbiased and thus gapless bilayer graphene. recent ex- 
perimental and theoretical studies 14, 15, l^, 
have revealed some interesting aspects of biased bilayer 
graphene. It has been shown that an electronic gap be- 
tween the valence and conduction bands opens up at the 
Dirac point and the low energy band acquires a Mex- 
ican hat dispersion relation by changing the density of 
charge carriers in the layers through the application of 
an external field or by chemical doping, which creates 
a potential difference between the layers. The presence 
of the potential bias transforms the bilayer graphene into 
the only known semiconductor with a tunable energy gap 
and may open a way for developing photodetectors and 
lasers tunable by the electric field effect. 

Under strong perpendicular magnetic field, experimen- 
tal results have shown that biased bilayer graphene ex- 
hibits a pronounced plateau at zero Hall conductivity 
<72:y=0, which is absent in the unbiased case and can 
only be understood as due to the opening of a sizable gap 
between the valence and conduction bands [l^. Tight- 
binding calculations have shown that the existence of 
such a gap can have a si gnif icant effect on the Landau 
level (LL) spectrum 



the presence of disorder in biased bilayer graphene is still 
absent and hence greatly needed. 

In this work, we carry out a numerical study of the 
QHE in biased bilayer graphene in the presence of disor- 
der based upon a tight-binding model. The Hall con- 
ductivity near the band center exhibits a sequence of 
plateaus at a^y — ve^jh where v is an integer, as in 
the conventional QHE systems. The = plateau is ro- 
bust with its width proportional to the strength of bias, 
which is consistent with the experimental observation. 
We further investigate the effect of random disorder on 
the QHE by calculating the Thouless number j20j. In- 
terestingly, at an intermediate disorder strength, the en- 
ergy gap around Ej = disappears, which destroys the 
V = plateau, and the system undergoes a transition to a 
metallic state. In the strong-disorder (or weak-magnetic- 
field) regime, the QHE plateaus around the band center 
can be destroyed due to the float-up of extended levels 
toward the band center. The v — ±2 plateaus are the 
most stable ones, which disappear last. Furthermore, we 
have also calculated the energy spectrum for an open 
system (cylindric geometry), and performed numerically 
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While disorder effect is 
known to be crucial in the conventional QHE systems, 
in-depth understanding of the properties of the QHE in 



a Laughlin's gauge experiment [21|, 1201 by adiabatically 
inserting flux quantum to directly probe the quantum 
transport near the sample edges. No conducting edge 
states are observed in the i/ = energy gap, suggesting 
an insulating state with divergent resistivity. 

The paper is organized as follows. In Sec. II, we in- 
troduce the model Hamiltonian and formulas for the cal- 
culation. In Sec. HI, numerical results based on exact 
diagonalization and transport calculations are presented. 
Sec. IV concludes with a summary. 
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II. THE TIGHT-BINDING MODEL OF BIASED 
BILAYER GRAPHENE 



HI. RESULTS AND DISCUSSION 



We consider a bilayer graphene sample consisting of 
two coupled hexagonal lattices including inequivalent 
sublattices A, B on the bottom layer and A, B on the top 
layer. The two layers are arranged in the AB (Bernal) 
stacking [H, , where B atoms are located directly be- 
low A atoms, and A atoms are the centers of the hexagons 
in the other layer. Here, the in-plane nearest-neighbor 
hopping integral between A and B atoms or between A 
and B atoms is denoted by ^ab = IXb ~ ^^''^ 
interlayer coupling, we take into account the largest hop- 
ping integral between a B atom and the nearest A atom 
7~g = 7i, and the smaller hopping integral between an 

A atom and three nearest B atoms 7^g^ — 73. The val- 
ues of these hopping integrals are taken to be 70 — 3.16 
eV, 71= 0.39 eV, and 73 — 0.315 eV, the same as in 
Ref. [H. 

We assume that each monolayer graphene has to- 
tally Li, zigzag chains with atomic sites on each 
chain [25|. The size of the sample will be denoted as 
X Ly X Lz, where = 2 is the number of 
graphene monolayers stacked along the z direction. In 
the presence of an applied magnetic field perpendicular 
to the plane of the biased bilayer graphene, the lattice 
model in real space can be written the following form [l3| : 

- 71 E ^"'^S^bC.J - 73 E + h.c. 

+ '^{wi + ei)cjcj -I- (wj -I- e2)c|ci, (1) 



where 



), cj(cjg) are creating operators on A and 
B sublattices in the bottom layer, and c|(c^~), c^(c^~) 

lA jB 

are creating operators on A and B sublattices in the 
top layer. The sum denotes the intralayer nearest- 

neighbor hopping in both layers, stands for inter- 

layer hopping between the B sublattice in the bottom 
layer and the A sublattice in the top layer, and 
stands for the interlayer hopping between the A sublat- 
tice in the bottom layer and the B sublattice in the top 
layer, as described above. For the biased system the two 
layers gain different electrostatic potentials, and the cor- 
responding energy difference is given by Ag = £2 — 61 



where ei = — ^A„, and 62 



Ag. For illustrative pur- 



pose, a relatively large asymmetric potential Ag = O.O570 
is assumed. Wi is a random disorder potential uniformly 
distributed in the interval Wi G [-W/2,W/2]^o. The 
magnetic flux per hexagon (/) — 



o 



^ IS propor- 



tional to the strength of the applied magnetic field B, 
where M is assumed to be an integer. 



The Hall conductivity axy can be calculated by using 
the Kubo formula through exact diagonalization of the 
system Hamiltonian [l3|. In Fig. 1, the Hall conductiv- 
ity o'xy near the band center is plotted as a function of 
electron Fermi energy Ef for a clean sample {W = 0) 
of size = 96 X 24 X 2 with magnetic flux (j) = ||, for 
biased and unbiased cases. Since the Hall conductivity 
is antisymmetric about zero energy, we show it mainly 
in the negative energy region. As we can see, in the un- 
biased case, the Hall conductivity exhibits a sequence of 
plateaus at axy = ve^ /h, where v = kgs with k an integer 
and gs = 2 due to double- valley degeneracy d, (the 
spin degeneracy will contribute an additional factor 2, 
which is omitted here). The transition from the 1/ — —2 
plateau to = 2 plateau is continuous without a = 
plateau appearing in between, so that a step of height 
Ae^/h occurs at the neutrality point. However, when a 
bias is applied, the valley degeneracy is lifted due to the 
different projection natures in the two layers of the LL 
states in the K and K' valleys. The valley asymmetry 
has a strong effect on the LLs near zero energy, where 
the charge imbalance is saturated. As a consequence, 
the Hall conductivity is quantized as axy — ve'^/h, where 
v — kgs with k an integer and gs = 1 for each LL due 
to the split of double- valley degeneracy [l^. With each 
additional LL being occupied, the total Hall conductivity 
is increased by e'^/h. Around the particle-hole symmet- 
ric point Ef — 0, a pronounced plateau with axy = 
is found, which can only be understood as due to the 
opening of sizable gap, Ag, between the valence and con- 
ductance bands. The emerged zero Hall plateau is ac- 
companied by a huge peak in the longitudinal resistivity 
Pxxi indicating an insulating state. This behavior has 
been observed experimentally 15|. It implies that a di- 
verging pxx at the particle- hole symmetric point Ej — 0, 
in striking contrast to all the other Hall plateaus, where 
Pxx vanishes as same as in ordinary QHE. 

Now we study the effect of random disorder on the 
QHE around the band center in the biased bilayer 
graphene based upon the calculation of the Thouless 
number. In Fig. 2, the Hall conductivity axy and Thou- 
less number g around the band center are shown as func- 
tions of Ef for three different disorder strengths and a 
relatively weak magnetic flux </» = |f ■ In Fig. 2a, the cal- 
culated axy and Thouless number g at a, weak-disorder 
strength W = 0.2 are plotted. Clearly, each valley in the 
Thouless number corresponds to a Hall plateau and each 
peak corresponds to a critical point between two neigh- 
boring Hall plateaus. We will call the central valley at 
Ef — the — valley, the first one just above (below) 
it the ly = —1 {ly = 1) valley, the second one the v = —2 
{ly — 2) valley, and so on, as same as the Hall plateaus. In 
Fig. 2b, the Hall conductivity axy and Thouless number 
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FIG. 1: Hall conductivity near the band center of unbiased 
and biased bilayer graphene with (f) = ^ . The disorder 
strength and sample size are set toW = and A*' — 96 x 24 x 2. 
Here, the spin degree of freedom has been omitted. 
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g for a relatively strong-disorder strength W = 0.6 are 
plotted. We see that the plateaus with ±2, ±6 and ±10 
remain well quantized, and the other plateaus become 
indiscernible, because of their relatively small plateau 
widths. With increasing W, higher valleys in the Thou- 
less number g (with larger are destroyed first, indi- 
cating the destruction of the corresponding higher Hall 
plateau states. When W — 2.0, all the plateaus except 
for the I' — ±2 ones are destroyed (see Fig. 2c). The 
last two plateaus v = ±2 eventually disappear around 
W ~ 3.2. Thus we observed that the destruction of the 
QHE states near the band center are due to the float-up 
of extended levels toward zero energy. 

In Fig. 3a, we show the Hall conductivity 
a function oi Ef for a relatively strong magnetic flux 

= and three different system sizes TV = 24 x 12 x 2, 
TV = 48 X 24 X 2, iV = 96 X 24 X 2 at disorder strength 
W = 2.0. We can see that at this disorder strength, the 
transition from v = —2 plateau to u = 2 plateau becomes 
continuous. With increasing the system size, the width 
of the plateau v = ±2 remains nearly unchanged. The 
region around the zero energy of Fig. 3a is enlarged in 
Fig. 3b. For comparison, we also show the results for 
the unbiased case, which clearly demonstrate the contin- 
uous behavior between the v = —2 plateau to the v = 2 
plateau in both cases. This behavior indicates a metal- 
lic state occurs around zero energy, which is essentially 
caused by the strong coupling between the two Dirac LLs 
due to disorder scattering. 

We now investigate the evolution of the edge states in 
an infinitesimal electric field by performing the Laugh- 



lin's gauge experiment 2l|, |22|. A periodic boundary 



condition in the y direction and an open boundary con- 
dition in the x direction are imposed to the system. The 



FIG. 2: Calculated Thouless number and Hall conductivity 
for = fl and three different disorder strengths, which are 
averaged over 400 disorder configurations. Here, the sample 
sizes are taken to be = 96 x 48 x 2 and A'' = 96 x 24 x 
2 in the calculations of the Thouless number and the Hall 
conductivity, respectively. 
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FIG. 3: Calculated Hall conductivity for biased and unbi- 
ased bilayer graphene with magnetic flux (f> = ^ at disorder 
strength W = 2.0 for three different system sizes. 



system can thus be considered as a cylinder. When the 
flux Oy{t) threading the cylinder is adiabatically turned 
on from Oy{0) = to dy{t)—2TT, which is equivalent to 
applying a weak electric field along the y direction 



Eyit) 



1 d9{t) 
Ly dt 
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. By diagonalizing the Hamiltonian Eq.(l) under the 
open boundary condition along the x-axis, at 200 dif- 
ferent 6y, the eigenenergies En of the system are ob- 
tained. Fig. 4a shows the calculated energy spectrum 
as a function of 9y for a clean sample {W = 0) at sys- 
tem size N = 96 X 24 X 2. Note that 9y ^ and 
9y=2'K are equivalent, as the system hamiltonian is peri- 
odic H(0y = O)=H(0j,=27r). We first examine the energy 
spectrum corresponding to the v = —2 QHE plateau. 
We observe that with changing 9y^ the energy levels in 
the plateau region cross each other, which correspond 
to two conducting edge channels in accordance with the 
quantized Hall effect. For example, we choose Fermi en- 
ergy Ef = O.I70. For 9y = 0, in the ground state, all 
the single particle states below Ef are occupied, whereas 
unoccupied above Ef. Upon insertion of the flux quan- 
tum, the two occupied states below Ef are pumped onto 
states above Ef indicated by the arrow, which causes two 
electrons transferred across from one edge to the other, 
corresponding to the quantized Hall conductivity with 
axy — 2e^//i, as shown in Fig. 4b. However, there are 
no such conducting edge states near Ef = 0.0, where 
the V — Q plateau is found. Clearly, a true spectrum 
gap shows up corresponding to a trivial insulating phase, 
which results in zero net charge transfer, and the current 
carried around the ribbon loop is zero. 

Now we consider the disorder effect. Fig. 5a shows the 
results for a randomly chosen disorder configuration for 
W = 2.0 at system size = 96 x 24 x 2. We can see that 
the energy gap around Ef — Q disappears. This behavior 
indicates that the transition from v = —2 plateau to 
u = 2 plateau becomes continuous, as shown in Fig. 5b. 
In contrast, if we choose an arbitrary Fermi energy in 
the V = ±2 plateau regions, e.g., Ef = O.I670, there are 
always two electrons transferred across from one edge to 
the other. Before the v = 2 plateau is destroyed by the 
disorder, the Ef ~ point becomes metallic. 



..(a) biased bilayer,.--"^^^^^^ 



(b) 



9 <',,(e-=/h) 

FIG. 4: (a) Energy levels of biased bilayer graphene with an 
open boundary in the x direction, as a function of the twisted 
boundary phase Oy in the y direction, (b) Hall conductivity 
near the band center for W = 0. Here (j) — ^ and = 
96 X 24 X 2. 
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IV. SUMMARY 

In summary, we have numerically investigated the 
QHE in biased bilayer graphene based on tight-binding 
model in the presence of disorder. The experimentally 
observed unconventional QHE is reproduced near the 
band center, where the Hall conductivity is quantized as 
o'xy = ve^ jh with v being any integer, including v = 0. 
The V = plateau around Ef = is due to the opening of 
sizable gap between the valence and conductance bands, 
which is absent in the unbiased case. By performing nu- 
merically a laughlin's gauge experiment, we have found 
that there are no conducting edge states in the ly — 
plateau region, in contrast to the ly ^ plateaus, where 
energy levels across each other, resulting in charge trans- 
fer between the edges and charge accumulation at the 
edges. However, at an intermediate disorder strength, the 



FIG. 5: (a)Energy levels of unbiased bilayer graphene as a 
function of twisted phase 9y. (b) (b) Hall conductivity near 
the band center for W = 0. Here (/> = || and = 96 x 24 x 2. 



energy gap around Ef — disappears, which indicates 
that the transition from 1/ = —2 plateau to v — 2 plateau 
becomes continuous, in agreement with the calculated re- 
sults of the Hall conductivity. Furthermore, we show that 
with increasing disorder strength, the Hall plateaus can 
be destroyed through the float-up of extended levels to- 
ward the band center and higher plateaus disappear first. 
At a strong-critical-disorder strength W = Wc — 3.2, the 
most stable QHE states with = ±2 eventually disap- 
pear, which indicates a transition of all the QHE phases 
into an insulating phase. 

Acknovifledgment: This work is supported by the US 
DOE grant DE-FG02-06ER46305 (LJZ, DNS), the NSF 



5 



grant DMR-0605696 (RM, DNS). We thank the KITP 
for partial support through the NSF grant PHY05-51164. 
We also thank the partial support from the State Schol- 
arship Fund from the China Scholarship Council, the Sci- 
entific Research Foundation of Graduate School of South- 
east University of China (RM) , the doctoral foundation of 
Chinese Universities under grant No. 20060286044(ML), 
the National Basic Research Program of China under 
grant Nos.: 2007CB925104 and 2009CB929504 (LS), and 
the NSF of China grant No.: 10874066 (LS). 



[1] K. S. Novoselov, E. McCann, S. V. Morozov, V. I. Falko, 

M. I. Katsnelson, U. Zeitler, D. Jiang, F. Schedin, and 

A. K. Geim, Nat. Phys. 2, 177 (2006). 
[2] R. V. Gorbachev, F. V. Tikhonenko, A. S. Mayorov, 

D.W. Horsell, and A. K. Savchenko, Phys. Rev. Lett. 

98, 176805 (2007). 
[3] S. V. Morozov, K. S. Novoselov, M. L Katsnelson, F. 

Schedin, D. C. Elias, J. A. Jaszczak, and A. K. Geim, 

Phys. Rev. Lett. 100, 016602 (2008). 
[4] E. A. Henriksen, Z. Jiang, L. C. Tung, M. E. Schwartz, 

M. Takita, Y. J. Wang, P. Kim, and H. L. Stormer, Phys. 

Rev. Lett. 100, 087403 (2008). 
[5] E. McCann and V. L Fal'ko, Phys. Rev. Lett. 96, 086805 

(2006). 

[6] J. Nilsson, A. H. Castro Neto, N. M. R. Peres, and F. 

Guinea, Phys. Rev. B 73, 214418 (2006). 
[7] J. G. Checkelsky, L. Li and N. P. Ong, Phys. Rev. Lett. 

100, 206801 (2008). 
[8] Y. Hasegawa and M. Kohmoto, Phys. Rev. B 74, 155415 

(2006). 

[9] D. A. Abanin, K. S. Novoselov, U. Zeitler, P. A. Lee, A. 



K. Geim and L. S. Levitov, Phys. Rev. Lett. 98, 196806 

(2007) . 

[10] E. V. Gorbar, V. P. Gusynin, V. A. Miransky and L A. 

Shovkovy, Phys. Rev. B 78, 085437 (2008). 
[11] H. Min and A.H. MacDonaId,Phys. Rev. B 77, 155416 

(2008) . 

[12] E. V. Castro, K. S. Novoselov, S. V. Morozov, N. 
M. R. Peres, J. M. B. Lopes dos Santos, J. Nils- 
son, F. Guinea, A. K. Geim, and A.H. Castro Neto, 
larXiv:0807. 3348 ^2008). 

[13] R. Ma, L. Sheng, R. Shen, M. Liu and D. N. Sheng, 
arXiv:0810. 1494^2008). 

[14] T. Ohta, A. Bostwick, T. Seyller, K. Horn and E. Roten- 
berg. Science 313, 951 (2006). 

[15] E. V. Castro, K. S. Novoselov, S. V. Morozov, N. M. R. 
Peres, J.M.B. Lopes dos Santos, J. Nilsson, F. Guinea, 
A. K. Geim, A. H. Castro Neto, Phys. Rev. Lett. 99, 
216802(2007). 

[16] J. B. Oostinga, H. B. Heersche, X. Liu, A. F. Morpurgo, 

and L. M. K. Vandersypen, Nat. Mater. 7, 151 (2008). 
[17] F. Guinea, A. H. Castro Neto, and N. M. R. Peres, Phys. 

Rev. B 73, 245426 (2006). 
[18] H. Min, B. Sahu, S. K. Banerjee, and A. H. MacDonald, 

Phys. Rev. B 75, 155115 (2007). 
[19] E. McCann, Phys. Rev. B 74, 161403(R) (2006). 
[20] J. T. Edwards and D. J. Thouless, J. Phys. C 5, 807 

(1972); D. J. Thouless, Phys. Rep. C 13, 93 (1974). 
[21] R. B. Laughlin, Phys. Rev. B 23, 5632 (1981). 
[22] B. L Halperin, Phys. Rev. B 25, 2185 (1982). 
[23] S. B. Trickey, F. Muller-Plathe, and G. H. F. Diercksen, 

Phys. Rev. B 45, 4460 (1992). 
[24] K. Yoshizawa, T. Kato, and T. Yamabe, J. Chem. Phys. 

105, 2099 (1996); T. Yumura and K. Yoshizawa, Chem. 

Phys. 279, 111 (2002). 
[25] D.N. Sheng, L. Sheng, and Z.Y. Weng, Phys. Rev. B 73, 

233406 (2006). 



